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Numerical solutions for: the integral curves of the velocity field (streamlines), the density con- 
tours, and the accretion rate of a steady-state flow of an ideal fluid with p = Kn'' equation of 
state orbiting in a core-dipole-shell system are presented. For 7 7^ 2, we found that the non-linear 
contribution appearing in the partial differential equation for the velocity potential has little effect 
in the form of the streamlines and density contour lines, but can be noticed in the density values. 
The study of several cases indicates that this appears to be the general situation. The accretion 
rate was found to increase when the constant 7 decreases. 
PACS numbers: 47.15. Hg, 04.25.Dm, 47.75. -f-f, 97.60. Lf 

I. INTRODUCTION 

In a recent work 0] we study the streamlines and density contours of a stationary fluid in the presence of either 
a rigid sphere or a black hole both with a distant shell of matter modeled by a dipolar field. To fix ideas see Fig. 1 
where this situation is depicted. The substitution of a halo or shell of matter by a fixed axisymmetric field is an old 
and common method in astrophysics, for example to simulate galaxy halos Numerical experiments confirm that 
this is usually a good approximation ||^. This method can be used in other physical situations where a core and a 
distant shell of matter are present. 

In we used for the fluid the stiff equation of state, i.e., a polytropic, p = Kn'' with 7 = 2. This stiff equation of 
state over simplifies the partial differential equation for the velocity potential field $ (it becomes linear) . In this case 
the sound velocity in the fiuid is equal to the speed of light. Therefore the stiff equation of state represents a limit 
situation not easily encountered in usual physical situations. Abrahams and Shapiro 0| studied the more realistic 
situation, a fiuid with polytropic equation of state with 1 < 7 < 2 in the presence of either a rigid sphere or a black 
hole. The sound velocity in this case is less than the speed of light and the partial differential equation for $ turns to 
be nonlinear. 

In the present work, following Abrahams and Shapiro |Q, we study a more realistic situation than the one examined 
in ||l|]. Now we compute the streamlines and baryon density contours for a fiuid with polytropic equation of state 
with 1 < 7 < 2 in the presence of either a rigid sphere or a black hole both with dipolar halos. In other words, now 
we consider the same situation as before, but for a general polytropic fiuid. We consider also the accretion rate of 
particles into the black hole at its dependence with 7. Furthermore, we complete the work of Abrahams and Shapiro 
Q by examining the black hole case with 1 < 7 < 2 and asymptotic constant velocity (case not studied in the quoted 
work). In particular, we study streamlines, density contours and accretion rate. 

In section II we present the basic equations that describe potential fiows and the nonlinear equation for the velocity 
potential for a polytropic equation of state of the form p = Kri^ . Section III is divided in three sub-sections. In 
sub-section III-A we introduce the metric that represents a central core with a distant dipolar shell of matter. This 
field is the first approximation to represent distant matter, like halos and rings. In sub-section III-B we present the 
numerical method used to solve the nonlinear partial differential equation presented in section II. In sub-section III-C 
we show some results and study the behavior of the streamlines and density contours of the fiuid in the presence of 
either a rigid sphere or a black hole both with a dipolar halo. Also we compute the accretion rate of particles falling 
into a black hole with and without dipolar halo. Finally, in section IV we summarize our results. 



II. BASIC EQUATIONS 

The solution for potential fiow of the equation of motion derived from the energy-momentum tensor for an ideal 
fluid, T^^ = {p + p)U^Uy +pgf_,^, is 
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hUa = (1) 

where h — {p + p)/n denotes the enthalpy per baryon, n the baryon density, p the pressure of the fluid, p = po + e 
the total energy density, po the rest mass energy density, e the internal energy density, represents the fluid four- 
velocity, and $ the scalar potential. Our conventions are G = c = f , metric with signature +2. Partial and covariant 
derivatives are denoted by commas and semicolons, respectively. 

We recall that, by definition, an ideal fluid is adiabatic i.e. the entropy is constant along the world line of each 
fluid element. If in an instant t the entropy is constant throughout the volume of the fluid (isentropic case) then, for 
all time and any subsequent motion of the fluid, it retains everywhere the same constant value. In that case we find 
the solution for potential flows Hence, only isentropic flows can have potential flows 

From Eq. (Q) and the baryon number density conservation equation, 

{nUn-,^ = 0, (2) 
we obtain the differential equation for the velocity potential, 

□ a>+[ln(^)],„ci>^" = 0, (3) 

where, — ['^/~gg^'^^^p]^^/^—g. This equation is in general nonlinear and depends on the fluid equation of state. 
If we consider a polytropic equation of state p = Kn^ and use the first law of thermodynamics for isentropic flows, 
d(p/n) +pd(l/n) — Tds, (s is the entropy per baryon and T the temperature), we find 

P = Po + 7- (4) 

7-1 

Now, assuming that the flow is relativistic we can neglect the rest mass energy density po with respect to the internal 
energy density Kn"^ /{'y— 1), i.e. the flow satisfies a barotropic equation of state, p — {'y— l)p- In this case, the baryon 
number density can be written in terms of the enthalpy as 

n — nhi-^ , (5) 

/ _,\ 1/(7-1) . 
where k = ( ^^^^^ j . Then, equation (^) can be written in the form 

□ $ + ^—^ [In h(<^>)] = 0, (6) 

7-1 

The simplest case is found when 7 = 2, Eq. reduces to a linear equation. In this case the barotropic sound 
speed, defined as dP/dp = cl, is equal to the speed of Hght, i.e. we have a stiff equation of state. In the general case 
(1 < 7 < 2), the differential equation is nonlinear and the barotropic sound speed is = (7 — 1) is less than the 
speed of light. 

The normalization condition, UaU" = — 1, gives us a relation between the enthalpy and the scalar field, 

h = (7) 

This relation will be also useful to determine the baryon number density. 

In this section we have followed the work of Moncrief |^] which is in accord with Tabensky and Taub [|| that 
consider a constant barotropic sound speed, = (7— 1). We note that Abrahams and Shapiro Q consider a variable 
polytropic sound speed. 

III. POTENTIAL FLOWS IN BLACK HOLES AND RIGID SPHERES WITH DIPOLAR HALOS 

A. The metric 

To incorporate a dipolar field in the Schwarzschild metric we consider the solution found in |^ where authors model 
the intermediate vacuum between a core and a distant shell of matter. This core-shell system was found solving the 



vacuum Einstein equations for a general static axially symmetric metric (Weyl solution |10|). One of these equations 
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can be solved in terms of Legendre polynomials increasing with the distance. These shell-like structures can be used 
to model many situations of interest in astrophysics, such as: supernovas, nebulae, and galaxy halos that have a core 
and an exterior shell of matter. As a first approximation we will consider only the first term in the expansion, the 
dipolar term that in Newtonian gravity corresponds to a constant force. By letting in the solution presented in |^ 
the quadrupole and octopole moments equal zero, we obtain 



dv 



2 1 



+ + 1)^(1 -i.2)e-2P™V 



(8) 
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+ (z + m)2 + p'^ + {z — m)2 
u > 1, 



p'^ + {z + m)2 — -^/p^ + (z — m)2 
-1 < -y < 1, 



(9) 



where (r, 0, (/?), (u,v,ip) and {p,z,ip) are spherical, prolate spherical, and cylindrical coordinates, respectively. The 
metric (||) represents a monopolar core in the presence of an external dipolar field {T>) that is associated to a distant 
shell or halo of matter. As we said before, this field can be seen as the first approximation to model an external 
concentration of matter like halos and rings. 



B. Numerical method 

To solve equation (^) in the space-time with metric (^) we assume: a) That fluid is stationary, i.e., the function 
$ depends on time only through the addition of —at, where a is a constant related to the zeroth component of the 
velocity, b) That due to the axial symmetry of the metric the potential $ does not depend on the variable (p, and c) 
That the fluid is a test fluid, i.e. the metric does not evolve and it is given a priori. Also we put m = 1. 

First we notice that due to the acceleration of the fluid produced by the dipole, the stationary condition leads to a 
region of instability away from the core (black hole or hard sphere). So we have a region of stability $7 that is depicted 
m Fig. g. The size of this region is very sensitive to the values of the dipolar field (I?) that accelerates the fluid. 
The stationary criteria is checked using the fact that baryon density ^ is a positive defined quantity |0| . With the 
above mentioned assumptions Eq. (|^) reduces to an elliptical differential equation with an inner boundary conditions 
near the black hole or rigid sphere and an external boundary condition (asymptotic condition). 

To solve Eq. (^) we use a computational code based on a marching method with a second order precision, five 
point, finite difference. The numerical grid is evenly spaced in r and in 9. The marching method is the Stabilized 
Error Vector Propagation (SEVP) which is very efficient in solving separable and non-separable elliptic equations. 
The main ingredient of the method is a clever superposition of a particular solution of the elliptic equation with and 
a homogeneous solution. For a detailed discussion of the method, see for instance [pl] , p^ . 

In the case of a rigid sphere, for the inner boundary condition, the usual condition of zero normal velocity in the 
surface of the sphere is employed. In general, the fluid velocity must be equal to the corresponding component of the 
velocity of the surface. Since usual stars have gaseous surfaces (not hard), this condition does not describe a typical 
flow around a star. In special astrophysical situations, like relativistic flows around a neutron star, this condition can 
be valid. For a black hole we use the condition that the fluid number density particle must remain finite on the black 
hole horizon ||l^. This leads to the numerical condition jj] that near u = 1 (the black hole horizon). 



dr* 



u-l 



= 0, (10) 



where r* = (u + 1) + 2 In {j^') is the tortoise radial coordinate [Q. With the assumption that (f>.f = —a < the 
condition ( p^ can also be satisfied in our case and it will be taken as the inner boundary condition. For numerical 
applications of this condition see [p],F5| . 
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For the outer boundary condition, with help of the geodesic equation for the metric (|^, we can find a characteristic 
value for the scalar field potential 

where the integration is performed along the lines of constant v, usually between the surface of the sphere (or black 
hole horizon) and 60 Schwarzschild radius. Later, with the results obtained for the scalar field $, we check the 
stationary condition to see the value of the radius in which the system becomes non-stationary. The constant k is the 
fluid particle energy function. 

To solve the non-linear equation (^) we first put the non-linear term equal to zero and calculate the linear part 
of the equation to find $. This solution is used as an initial guess for the non- linear problem. Then, from (^ we 
compute the fluid enthalpy. Finally, with this information we compute the non-linear term that is introduced in the 
non-linear equation to find new values of The process is repeated until the sum of the fractional change in the 
enthalpy for all interior points of the grid in one iteration is less than a certain error e. Usually less than 10 iterations 
were required to reach an error of e < 10^^. Obviously, for the linear case we need only one iteration. The number of 
iterations depends on the value of 7. For 7 « 1 the nonlinear factor tends to infinity and the program diverge. When 
7 = 1 the fluid is pressureless (dust) and the fluid flow is geodesic, i.e., no longer obeys (^). 

The code was tested for the case j — 2 using the analytical results for the steady flow of a fluid in the presence 
of either a hard sphere jl^ or a black hole ||l^ . In this case the numerical solution agreed with the exact solution 
within an error better than 1.5% for the radial velocity and 1% for the angular velocity. Since we are interested in 
qualitative aspects of the fluid behavior this accuracy is sufficient for our purposes. 

Another aspect of the fluid dynamics is the accretion rate of matter into a black hole. This accretion rate can be 
computed from (|^) and (||), we find 

N = - I nW^f^dS^ = - [ Kh^<^,ug''"^/^dudv, (12) 
Js Js 

where the integration is performed in a two-surface sphere centered on the black hole. The exact form of (^2|), for 
7 = 2, is known ||l^. It is TV = 167rM^nooC^^, where Uoo and are the asymptotic density and zeroth component 
of the 4-velocity, respectively. M is the black hole mass. We test our code to compute the accretion rate with this 
exact solution. We found an error less than 1%. 



C. Numerical results 

In Figs. ^ and ^, for the case of a rigid sphere with a dipolar halo, we show the numerical results for the streamlines 
and density contour lines of the baryon number density. The constants are: radius of the sphere 2.5 (Schwarzschild 
radius 2), dipolar strength V = 0.001, 7 = 1.75 and k ^ a = 1.25, the outer boundary condition was set to 60 
Schwarzschild radius. We see that the streamlines do not differ much from the linear case 0] and that the contour 
lines have the same qualitative behavior as in the linear case, (7 = 2). The numerical values for the baryon density 
are not the same due to the exponent 1/(7—1) presented in the relation between the enthalpy and the baryon density. 
To ensure a physical result, in solving (g), we require that enthalpy be positive in every iteration of our code. The 
baryon number density is a positive quantity. We can have a negative enthalpy solution of n = Kh^/^''~^\ but in this 
case we have a negative baryon density. 

For the same values of the constants, in Figs. || and || we show the numerical results for the streamlines and density 
contours for the case of a black hole with a dipolar halo. Again we see little difference compared to the linear case. 
These results tell us that the contribution of the non-linear term does not affect in a significant way the qualitative 
behavior of the fluid. 

For the case of a black hole without dipolc {V — 0) we also found that, for a constant asymptotic velocity, that the 
streamlines and density contours do not distinguish between the linear and non linear case. The scaled accretion rate 
for the black hole {N / k) increases when 7 decreases, e.g., for 7 = 1.9, 1.8, 1.7, the scaled accretion rate increases in 
6%, 13%, 24%, respectively, the value of the asymptotic velocity used in this case is Voa = 0.6, the same as in 

For 7 = 2, the linear case, we find that the value of the scaled accretion rate with or without dipolar halo differs by 
1%. We also find in this case, for the same values of Fig. ^, that the scaled accretion rate increases when 7 decreases, 
e.g., for 7 = 1.9,1.8,1.7, the scaled accretion rate increases in 8%, 18%, 33%, respectively. In the dipole case the 



4 



growth of the scaled accretion rate when 7 decreases is greater than in the case without dipole. In both cases the 
growth of the scaled accretion rate is due to the factor involving the enthalpy [cf. Eq. . 

As mentioned before, for an accelerated fluid it is difficult to have stationary flow. For the cases studied along the 
paper we found that for T) < 0.01 we can have a region of reasonable size (a ball greater than 45 Schwarzschild radius) 
where we have laminar flow. For D = 0.03 the radius of this ball shrinks to 5, and for V — 0.05 stationary flow does 
not exist. In summary, we believe that our results are representative of a generic situation when V < 0.01. 

Finally, we want to point out that the variation of $ in the systems black hole and rigid sphere with or without 
dipole for the linear and nonlinear cases is, in almost every interior point of the grid, less than 5%. 



IV. CONCLUSIONS 



We have solved numerically the streamlines and the density contours for the baryon number density of an ideal 
fluid in the presence of either a rigid sphere or a black hole both with dipolar halo. The fluid has barotropic equation 
of state, p = (7 — 1) with 1 < 7 < 2, that represents a more realistic physical situation than the linear case studied in 
. We see that the nonlinear term does not affect qualitative the form of the streamlines and density contours when 
compared to the linear case, this is due to the low variation of the scalar field between the linear and nonlinear cases. 
When compared with other situations studied in the literature this fact appears to be generic. Another interesting 
result is that the accretion rate of the cases studied increases when the constant 7 decreases. In the presence of the 
dipolar field this growth is bigger compared with the case without dipole. 

The cases of a rotating black hole and rigid sphere with halos modeled by multipole moments beyond the dipole is 
under active consideration by the authors. 
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FIG. 2. Numerical results for the fluid streamlines with 7 — 1.75 when an external dipolar field with D — 0.001 is present 
and a rigid sphere of radius r = 2.5 (Schwarzschild radius equal 2) is placed as an obstacle. We set the values of the constant 
A; = a = 1.25 . The axes are defined as X = rsm9 and Z = r cos 6, with r = u + 1 and cos 6 = v. 
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FIG. 3. Numerical results for the density contours {n/n) of the baryon number density when an external dipolar field with 
T) — 0.001 is present and a rigid sphere of radius r — 2.5 (Schwarzschild radius equal 2) is placed as an obstacle. The axes and 
the constants are defined as in Fig. ^. 




FIG. 4. Numerical results for the streamlines when an external dipolar field of value V = 0.001 is present and a black hole 
is placed as an obstacle. The black hole has radius r — 2 (Schwarzschild radius equal 2). The constants and the meaning of 
axes are the same of Fig. ^ 
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